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ABSTRACT 

Purpose: The purpose of the study is the quantitative assessment of Electrical Impedance Endotomography (EIE) 
for the specification of hardware systems. EIE is a modality of Electrical Impedance Tomography (EIT) where the 
electrodes are located on a probe placed in the middle of the region of interest. The absence of material boundary to the 
explored volume and the decrease in sensitivity away from the probe requires specific study. 

Material and methods: The method is the derivation of the equation linking explored medium's conductivity, the 
sensitivity distribution of the electrode patterns used for data collection and measuring system's noise and bandwidth. 
The assessment of EIE was achieved by means of simulations based on realistic data of conductivity and noise level. 

Results: The derived equation enabled the estimation of the current needed under realistic operating conditions 
corresponding to prostate imaging. The generalisation to other organs is straightforward. The image reconstructed from 
the simulated data and from bench experiments were in agreement and showed that the two selected drive patterns, fan3 
and adjacent, gave images of similar quality in absence of noise and that adjacent drive requires significantly higher 
measurement current. 

Conclusion: The study confirmed the feasibility of EIE with achievable hardware specifications. The derived 
equation enabled the determination of design parameters for the specification of hardware systems corresponding to any 
given application. The study also showed that EIE is more appropriate for tissue characterisation than for high speed 
imaging. © 2006 Biomedical Imaging and Intervention Journal. All rights reserved. 
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INTRODUCTION 

Electrical Impedance Tomography (EIT) produces 
images of a body from impedance data collected using 
surface electrodes. The advantages of this method in 
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medical applications include harmlessness, ease-of-use, 
high time resolution and specificity for tissue 
characterisation [1-4]. The direct problem obeys the 
second order differential equation V (oVu) = 0 with 
Neumann/Dirichlet mixed boundary conditions. The 
associated inverse problem is ill-posed [5]. The small 
number of measurements limits the spatial resolution of 
the reconstructed images. The reduction of sensitivity 
and resolution from periphery towards the centre [6] 
makes EIT not well suited for the imaging of small and 
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Figure 1 Sketch of the tip of a 16-electrode probe for EIE. The insulating core may consist of a cylindrical tube 
enabling to pass electrode leads. 



deep-located organs. The nature of the governing 
equation and the low number of possible measurements 
with a workable number of electrodes prevents 
significant improvement of spatial resolution. 

In EIE, the electrodes are located on a probe placed 
in the middle of the region of interest for local 
measurements. An EIE probe consists of multiple linear 
electrodes regularly spaced on the outer surface of an 
insulating core (Figure 1). EIE was developed for 
prostate imaging aimed at the evaluation of cancer 
treatment by therapeutic ultrasound [7-9]. In the 
proposed application, the impedance method is expected 
to compensate for the lack of specificity of ultrasonic 
imaging in cancer detection and the low acoustic contrast 
observed in the prostate between normal tissue and tissue 
treated by ultrasound [10-12]. This approach has been 
supported by the known significant conductivity 
differences between cancerous and normal tissue 
observed in various organs [13-18], including human 
prostate [19,20], and by recent studies reporting 
significant conductivity changes in tissue exposed to 
ultrasound energy [21,22]. More generally, EIE can 
potentially address a range of interstitial and intracanular 
measurements such as in oesophageal and vascular 
studies. 

The method, however, does not aim to compete with 
the radiological methods of prostate imaging. The 
objective is to derive a complementary technique for use 
in conjunction with ultrasound techniques. The objective 
is to exploit the specificity of impedance measurements 
to improve the characterisation of tissue before and after 
treatments with therapeutic ultrasound. 

EIE obeys the same governing equation as EIT but 
with different boundary conditions. The boundary profile 
is unknown and variable in EIT according to the 
morphology of the examined region and inter-patient 
variability while the surface bearing the electrodes is 
known by construction in EIE. It is obvious that the 
volume actually sensed by the probe is finite although 
there is no material boundary around the probe. The limit 



is the distance beyond which noise overrides the 
contributions of distant points. 

The determination of the computational domain is a 
key problem in EIE. The absence of tangible limit for the 
domain sensed by the EIE probe is one major difference 
compared with EIT. Therefore, the notions formerly 
investigated in EIT are still relevant in EIE, due to the 
same governing equation, but need to be revisited. For 
quantitative studies, the geometry of an EIE probe 
suggested the use of a model with axial symmetry. In this 
2D model, an infinitely long cylinder represents the core 
of the probe and infinitely long lines regularly spaced at 
the outer surface of this cylinder represent the electrodes. 
This model enabled the derivation of analytical equations 
for current density, field and potential created by EIE 
electrodes in a medium of homogeneous conductivity. 

According to the model, the field of a single 
electrode is proportional to 1/d, where d is the distance to 
the centre of the probe. Hence, the field of a pair of 
electrodes tends to vary as 1/d 2 for large distances to the 
probe. The consequence is a rapid fall-off in sensitivity 
of the order of 1/d 4 . The measurements in a semi-infinite 
medium present a certain similarity with the case of the 
rosette array of surface electrodes used in EIT for 
monitoring gastric function [23]. In both cases, the 
electrodes are grouped together, do not encircle the 
region of interest, and explore a semi-infinite medium. 
The measurements carried out with the rosette support 
the feasibility of EIE measurements. 

It was found that for a total number of 16 electrodes, 
the possible 4-electrode patterns using an adjacent pair of 
voltage electrodes can be sorted into 49 basic patterns 
from which any pattern can be obtained by symmetry 
and rotation. The association of the model with the lead 
field theory enabled the calculations of sensitivity maps 
for the basic 49 patterns. The study of these maps 
showed that the extension of sensitivity increases with 
the angular spacing of source electrodes. These maps are 
shown in the form of an animation (Figure 1 1). The drive 
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Figure 2 Reconstruction mesh used in this study. The central area corresponds to the insulating core of the probe. 

There are 14 layers and 64 angular sectors forming 896 trapezoidal pixels of constant profile but varying 
dimensions. 



pattern giving the largest sensitivity range was selected 
based on these sensitivity maps. 

The purpose of this study is the quantitative 
assessment of EIE as a whole including medium, 
electrodes and instrumentation. This differs from the 
previous studies, which were limited to the comparison 
of drive patterns to determine the widest sensitivity range. 
The novelty of the present study is to encompass all the 
components involved in EIE data collection. This was 
achieved by the derivation of an equation linking 
measurement, noise, magnitude of injected current, 
electrode sensitivity distribution, medium conductivity 
and conductivity contrast to be observed. The inclusion 
of noise enabled the calculation of the volume actually 
sensed by the probe. The study is supported by the 
comparison of adjacent and fan3 drive patterns using 
calculated data and images reconstructed from computer 
and experimental data. Although particular attention was 
given to prostate imaging, the study has been intended to 
enable generalisation to other applications and the design 
of hardware systems. 

DEFINITIONS 

Electrode patterns 

In this study, the measurements were carried out 
according to the 4-electrode technique with bipolar 
current patterns and differential voltage sensing with all 
four electrodes located on the probe. The sensing pair 
always consists of adjacent electrodes for hardware 
reasons including reduction of common mode signal. 
The number of patterns with adjacent voltage electrodes 
is N T = (N E -3) (N E -2)N E /2, where N E is the number of 



electrodes on the probe. The maximum number of 
linearly independent patterns is (N E -3)N E /2, as this was 
formerly demonstrated in EIT. In this study, it was found 
convenient to consider that the set of 4-electrode patterns 
used for data acquisition consists of the N E angular 
duplications of a basic set consisting of patterns 
comprising a given pair of voltage electrodes (pair 
arbitrarily denoted {0,1} in this study) associated to 
different pairs of current electrodes. 

Reconstruction mesh 

The reconstruction mesh consisted of N L concentric 
layers of N A trapezoidal pixels. The outer radius of the 
mesh was denoted R max . The vertices of a pixel were 
located on two circles of radii r n _i and r n , with r n _. < r n , 
r 0 =l and r NL = R max . The number of layers was N c = 14, 
the number of angular sectors was N A = 64. Hence, the 
number of pixels was N pix = 896. The mesh was designed, 
so that pixel dimensions were proportional to the 
distance from the origin (Figure 2). This was achieved in 
the setting of the radial increment equal to the length of 
the circular arc passing by the centre of the pixel. This 
condition can be written under the form of (1): 

rn=r n - l {*A+*WA-*) « 
In this mesh design, the resolution is better for 
pixels of higher sensitivity while the increase in pixel 
size with distance tends to compensate for the sensitivity 
decrease. The radius of the reconstructed domain was 
chosen according to the size of the domain explored by 
the probe. If the reconstruction domain is too small, 
significant elements would be ignored and attributed 
erroneously to pixels located inside the mesh. If the 
reconstruction domain is too wide, it would encompass 
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Figure 3 Construction of fan3 pattern. The current injection circuit, (S), is successively connected to the pairs of 
source electrodes. Fan3 comprises two groups of patterns. The first group uses electrode #3 successively 
associated to electrodes #9 to #15 (7 patterns). The second group (not shown in this figure for clarity) 
consists of symmetrical patterns versus the axis of symmetry of the voltage electrodes. In this group 
electrode #14 is successively associated to electrodes #2 to #8. One pattern {3, 14} is its own image in 
the symmetry, so that fan3 consists of 13 measurement patterns. 



points with negligible contributions. In this study, the 
reconstruction domain was determined by considering 
the diameters of urethral probes used in urological 
practice and the size of the prostate. This 3 -cm high 
organ is approximately conical in shape. It presents a 
base, an anterior, a posterior and two lateral surfaces. 
The base applied to the inferior surface of the bladder 
and the apex is directed downwards. The prostate is 
about 4x2 cm 2 at the base (2 cm in antero-posterior 
diameter). From these dimensions, the value of R max 
chosen was equal to 4 times the radius of the probe. This 
justified the use of 14-layer mesh. 

Sensitivity 

Sensitivity is a general concept to quantify the 
change in the measured signal when the conductivity 
within a given element, Ax, changes by Ao. The 
computation of sensitivity normally requires the 
resolution of the governing equation knowing the 
original and the perturbed distribution of conductivity in 
the medium. For small perturbations, the lead field 
theory yields a linear approximation that has been widely 
used in impedance imaging. This theory has enabled the 
derivation of a general expression of change, AZ X , in the 
measured impedance, Z x , due to the conductivity change 
in a given element Ax [24]. In the present study, the 
conductivity change was assumed to be uniform within 



the volume element. Hence, the general expression of 
AZ X , transforms into (2): 

AZ X «-Aa j Evo " E ' curr dz (2) 

(Ar) ^ volt I curr 

Evoit/Ivoit is the field of the voltage electrode in the 
initial medium. E' cun /I curr is the lead field of the current 
electrodes in the medium perturbed by Ao in element Ax. 
If the element volume and the conductivity change are 
small enough, it is convenient to consider that E' curr is 
approximately equal to the lead field of the current 
electrodes in the non-perturbed medium. Using Ohm's 
law, the voltage changes due to the conductivity change 
is given by (3), where I s denotes the measurement of the 
current injected across the source electrodes and o the 
initial conductivity value within the pixel: 

du = -AZ x xI s » (l s /(7)x(A<x/er)A r ^ 
with A f = -f J """ J ™ rr dt 

Al J vol, hurr 

Sensitivity A T is expressed in m" 1 in general (3D) 
and dimensionless in 2D (translationally uniform model). 
In the latter case, infinitely long lines model the 
electrodes and the measurement current in (2) is in Aim, 
so that the dimensions in (2) and (3) remain consistent. 

The sensitivity of all pixels and all electrode patterns 
used for data acquisition form the sensitivity matrix 
formed by N meas rows and N A xN L columns. The so-called 
"fan3" electrode pattern has a larger sensitivity domain 
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Figure 4 Two circuits for voltage-to-current conversion commonly used in impedance measurement. In both 
cases, the output current, IS, is equal to GVE/RS. In the Howland circuit (a), the voltage gain is 
G=R2/R1=R4/R3 and the noise gain is G+l. Using a difference amplifier (b), the voltage gain and noise 
gain are both equal to unity. 



than the other bipolar drive patterns tested: adjacent, 
diametric and fan4 [25]. 

The adjacent drive pattern has widely been used in 
EIT. It consists of 4-electrode patterns where both 
voltage and current electrode pairs consist of adjacent 
electrodes. Fan3 consists of 4-electrode patterns where 
the voltage electrodes are adjacent and source electrodes 
are of variable spacing (Figure 3). In fan3, the two 
source electrodes are separated by the symmetry axis of 
the voltage electrodes. The sixteen angular replications 
of these 13 patterns yields N meas = 208 measurement 
patterns. Similar definition applies to fan4 to fan8 
patterns. However fan3 was found to give the largest 
sensitivity range, so that the other fanX patterns were 
ignored in this study. 

NOISE 

Noise condition 

This section describes the sources of noise and 
presents the derivation of an equation for the 
measurement current. In the following, the term "signal" 
denotes the change in the voltage difference across the 
sensing electrodes in measuring the perturbed medium 
and the initial medium (5u in (3)). 

The correct measurement of a conductivity change, 
Ac, in an element implies that the sensitivity of this 



element is above the noise level. It was assumed in this 
study that the minimum conductivity change to be 
measured, denoted. | Act/ct | m j„, was the same for all pixels 
of the mesh. In all drive patterns, certain pixels have low 
sensitivity values due to either their distance to the probe 
or the local orthogonality of the lead fields. The 
contributions of such elements remain under the noise 
level for any realistic value of the measurement current. 
Hence, the noise condition used in this study was the 
following: any pixel of the mesh is sensed by at least one 
of the N E angular replications of any basic electrode 
pattern. This condition is really a minimal condition, for, 
if it were not satisfied, certain measurements would 
ignore certain pixels. This condition means that for any 
basic pattern, the contributions of at least N A /N E pixels 
per layer are above noise level. Hence, with a 16- 
electrode probe and the 64 angular sectors, the relevant 
parameter is then the fourth largest sensitivity magnitude 
in each layer of the mesh. 

Noise sources 

Three types of noise have been considered in the 
present study: electrode noise, electronic noise and 
current noise. Electrode noise's origin is electrochemical. 
In the frequency range of impedance measurements 
(f>l kHz), the general expression for such noise is 
similar to that of thermal noise [26-27]: 



JJossinet et al. Biomed Imaging Intern J 2006; 2(2):e24 



6 

This page number is not 
lor citation purpose 



» adja ->-fan3 o fan4 



8 




0.0 1.0 2.0 3.0 4.0 



Figure 5 Calculated measurement current (mA rms) needed to satisfy the noise condition of (9) with the mesh of 
Figure 2 under the conditions of the numerical example, section "Numerical application and |Ao/o| = 0.1 
(Section 4). The horizontal axis is the normalised distance (R pr0 bc = 1) from the origin to the centre of the 
pixel. 



e a = J4k¥m = Nf d Vai" Vb (4) 
with Nf el = 0.123nF(7nHz) ~' at 20°C 

where k is Boltzmann's constant, T absolute 
temperature (°K) and 9? the real part of the interface 
impedance and B system's bandwidth (Hz). 

For the quantification of the noise generated in the 
current injection circuit, the injecting circuit was 
considered as a voltage-to-current converter of 
transadmittance G/R s (S/m), where G is the voltage gain 
of the circuit and R s the resistance across which the 
feedback voltage is measured. Figure 3 shows two 
circuits commonly used for voltage-to-current 
conversion: Howland's circuit and differencing amplifier. 

The output current is Is=VexG/Rs. G denotes the 
closed loop voltage gain of the amplifier. Depending on 
the configuration of the circuit, the noise gain, G noise , can 
be different from G. Let be the noise figure (V/VHz) 
at the input of the circuit. The injected noise current, of 
rms value denoted i no i se , is then given by (5) where B is 
the bandwidth (Hz) of the system: 

1 noise = G noise N f, i4B I R s (5) 

This noise current produces an error voltage across 
the measured impedance Z x . Considering that the 
measured impedance is the quotient of the geometry 
factor g x to the mean conductivity of the medium, o m , the 
noise voltage due to the current source is given by (6): 

SivG^Nfjlg^jB/^Rs (6) 

The value of g x can be either calculated for each 
element and electrode pattern used by means of an 
appropriate model or derived from experimental 
measurements of impedance Z x , and medium's mean 
conductivity o m . Table 1 shows the values of g x 
calculated using the described 2D model and measured 



in vitro. The difference between measured and calculated 
values has been attributed to the dispersion of current 
streamlines at the extremities of finite electrodes [9] of 
length equal to the diameter of the probe as described in 
section "Experimental setup". 

The rms value of the amplifier input related noise 
voltage is given by (7): 

s v =Nf v jB (7) 

The noise figure Nf v is given in technical data 
sheets of operational amplifiers. Finally, using (4), (6) 
and (7), the total noise rms voltage superimposed to the 
signal Av x is given by: 

e T = 4l{s 2 d + sj +4]' 2 Jb (8) 
The coefficient a/2 accounts for the fact that the 
"signal" is the difference between two measurements. 

Numerical application 

The noise figure of the voltage amplifier is equal to 
13nV/VHz, which corresponds to standard op-amps 
usable in EIE. The contact impedance of one electrode of 
this probe in tap water of conductivity 0.039 S/m was 
400 Q. at the used measurement frequency of 8 kHz [8]. 
In a urethral probe, the electrode surface would be 
reduced by a factor of about 25 or less with respect to the 
mock-up probe used in bench experiments. The interface 
impedance would then become about 10 k Q in tap water. 
However, tissue conductivity is in general higher than 
that of tap water, so that the interface impedance will 
presumably be lower in tissue than in tap water. In the 
absence of literature data for urethral wall, data for blood 
vessel and prostate were considered instead. From the 
data published on the website of the Institute for Applied 
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Figure 6 Present version of the front end. The largest board comprises 16 electrode buffers, 16 differential 
voltage amplifiers. The other boards are two multiplexers boards for the selection source and sensing 
electrode pairs, the current source with the current measuring circuits and the control logic. The 
electrodes are connected to a 16-line linking the four stacked printed boards. Stacking boards enables 
flexibility for successive versions of the hardware system. 



Physics "Nello Carrara" [28], the magnitude of the 
calculated prostate admittivity ,(a*= o 0 +jcoee 0 ) varies 
from about 0.43 S/m to 0.58 S/m in the prostate and from 
0.28 to 0.33 for blood vessel in the range 
10 kHz - 1 MHz. Furthermore, the possible presence of 
urine, wetting urethra wall, would presumably tend to 
decrease the contact impedance. 

The presence of urine of higher conductivity than 
the surrounding tissue, can potentially affect the 
measurements. Besides the reduction of electrode contact 
impedance due to the wet urethral wall, the presence of 
urine could also cause adjacent electrodes to short circuit. 
It may be expected that the amount of urine present 
during the measurements would be limited by the 
preliminary draining of the urethra and the temporary 
obstruction of the lumen by the tip of probe. The shorting 
impedance would depend on the thickness of the 
conductive layer forming between the probe surface and 
the urethral wall. A possible protection measure would 
be to give electrode edges a slightly salient profile to 
locally increase the pressure to constrict or divide the 
conductive layer. This issue can only be solved by a 
practical measurement in situ. 

Conductivity values for urine range from 2.5 S/m to 
4.5 S/m [29, 30]. In this study, in the worst case the real 
part of the interface impedance of two electrodes in 
series was finally maximised by 2500 Q. The 
corresponding electrode noise figure, Nf d , is then about 
6 nVVHz. The noise figure due to the current injecting 
circuit can be derived from (6). With G n0 ise = 2, 
Nfi = 13nV/VHz, |g x | = 0.333 (maximal value in Table 1), 
R s = 1 kQ and o m = 1 S/m, the noise figure due to the 
current source is about 0.009nV/VHz, which is negligible 
compared with the other sources of noise. The total noise 



figure in this numerical example is then about 
14.3nV/VHz. 

Noise condition and measurement of current equation 

The condition for the correct measurement of the 
contribution of an element is that the contribution, 5u, of 
this element is above noise level. Assuming that the 
distribution of noise amplitude is Gaussian, one may take 
3.09xNf as arbitrary noise threshold, with the risk of 
0.002 for the instantaneous noise voltage, which is 
outside the interval ±3.09xN f . Using (4), (5), (6) in (8) 
and grouping the terms corresponding to the probe, the 
instrumentation and the medium finally gives the 
expression of the rms value of the injected current 
satisfying the above noise condition: 

/ s >(3.09V2/|A x |)x 

[[Nf d Mj + {G noise Nfig, l° m Rs) 2 + Wv J 2 x 

(9) 

The mean conductivity of the medium, o m , depends 
only on the explored medium, the conductivity change 
a/a on the observed phenomena (tumour, treatment). 

DESIGN PARAMETERS 

The purpose of this section is to determine the 
realistic design values for incorporation into (9) of the 
mean conductivity, c m , and minimal conductivity change 
[Ao/p[min- The mean conductivity value determines the 
magnitude of the measured impedance Z x , and 
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Figure 7 Singular values of the sensitivity matrices for diametric, adjacent, fan3 and fan4 bipolar drive patterns. 
The rank was 76 for diametric drive and 104 for all other drive patterns. 



consequently that of a pixel's contribution, 8u. For the 
prostate, Dawson [31] reports a value of 0.4 S/m for 
studies at 60 Hz. This value is close to the static 
conductivity o 0 of the prostate given by the IF AC 
internet resource [28]. The conductivity values calculated 
in section "Numerical application" from this resource for 
normal prostate tissue (0.43 to 0.58) suggested 0.5 S/m 
as design value of o m in (9). 

The design value of |Ao/o| m j n corresponds to the 
smaller change due to either the presence of cancer tissue 
or the treatment of a tissue by therapeutic ultrasound. 
Several sources of data enable the estimation of Ao/o in 
presence of cancerous tissue. Blad [32] has proposed a 
general conductivity ratio between normal tissue and 
cancer tissue of about 0.69 (|Ao/o| =0.31). Conductivity 
ratios of 0.72 and 0.78 at 16 kHz and 125 kHz, 
respectively, were observed between carcinoma and 
glandular tissue in excised breast tissue samples [17, 18]. 

Dunning tumour in a Copenhagen rat is a commonly 
used model for human prostate cancer [33]. The 
admittance of growing AT2 Dunning tumours was 
monitored during 21 days [34]. The conductivity was 
estimated by modelling the tumour with a cylindrical 
segment of length equal to its diameter with four equally 
spaced electrodes. This yielded a rough estimate of 
tumour conductivity of about 0.2 S/m at 9 kHz and 
0.4 S/m at 1 MHz. Using the figures and the values for 
reference prostate tissue of section "Numerical 
application", the coarse estimates of conductivity ratios 
are about 0.47 and 0.69 at the two considered frequencies. 

Lee [19] carried out impedance measurements at 
100 kHz, 1 MHz, 2 MHz and 4 MHz in prostates ex-vivo 
using bio-impedance needles. The conductivity was 
smaller in cancer tissue than in prostate tissue at 100 kHz, 



1 MHz and 2 MHz with conductivity ratios of 0.86, 0.92 
at and 0.8, respectively. 

Smith, by measuring eddy currents using a magnetic 
coil at 2.14 MHz, compared the conductivity values in 
Dunning tumours from G, AT2 and AT3 lines [20]. The 
conductivity was lower in AT2 and AT3 tumours (0.22 
and 0.24 S/m) and G line (0.33 S/m) than in control 
tissue (0.35 S/m). From the above data, the value of 10% 
was taken as representing the minimal change |Ao/o| m i n 
resulting from the presence of prostate cancer. 

The energy deposited in tissue by therapeutic 
ultrasound produces the irreversible necrosis of the tissue. 
In vitro experiments showed noticeable changes in a 
tissue's impedance. Changes larger than 20% were 
observed in tissue samples exposed in vitro to high 
energy ultrasound [21] and muscle tissue samples [22]. 
As these values were larger than the value for cancer 
tissue, the latter (10%) was finally taken as design value 
of |Ao/o| min for incorporation into (9). Figure 4 shows the 
plot of the measurement current satisfying the noise 
condition (9) with the value of section "Numerical 
application" and |Ao/o| min = 0.1. In this plot, the 
parameter varying with distance is A x , the fourth largest 
sensitivity value calculated for each layer of the mesh 
described in Figure 2. 



SIMULATION 

This section describes the simulation of an EIE 
application using calculated and experimental data. The 
conductivity ratio (o/o 0 ) was set to 1 . 1 for the simulation 
of biological conditions and to zero to simulate the 
plastic rods used in bench experiments. The model 
described below was used for the simulation of two 
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Figure 8 Images reconstructed from calculated noiseless data for adjacent and fan3 drive patterns. The simulated 
conductivity perturbations (Ag/g=0. 1 , radius r = 0.2) were centred at d = 2 (top row) and d = 3 (bottom 
row). Regularisation coefficient was X = I0' n for all four images. 



conductivity perturbations, of radius 0.2 and 0.3 and 
centred on the Oy axis at (0,2) and (0,3), respectively. 
The experimental data were collected using the bench 
model described in section "Experimental setup". 

Software model 

The signal, the change in the measured potential 
difference across a pair of sensing electrodes, was 
calculated using the 2D software model developed for 
the project [35]. In this model, infinitely long lines 
represent the electrodes and all quantities are assumed 
constant by translation along one direction. This models 
yields analytical equations for electric field and potential. 
The conductivity perturbations were assumed to be 
infinitely long cylinders parallel to the axis of the probe 
and projecting on the calculation plane as circular disks 
of conductivity g=Go+Ao. The voltage change at the 
sensing electrodes was calculated using the image theory 
[36] considering the series of images of the initial source 
electrodes in the perturbing cylinder and in the probe. 
This forms sequences of sources with rapid convergence 
of potential and electric field. 

The addition of noise to calculated data was 
achieved according to (1 1): 

^A k ) = ( I meas/ cr m) >< 

(M calc (k) + N L x\\M calc \\ 2 xX(k)) 



1 1 A.u 1 1 2 is the 2-norm of the data vector calculated for 
Is/Om equal to unity. N L denotes the noise level 
(dimensionless) and X a normal Gaussian variable. 
Hence, assuming that the variance of sample X k is equal 
to the variance of the distribution, the relation between 
total noise, e t and added noise is given by (12): 

N L=s T a m l{l meas \/S Mcalc \ 2 ) (12) 

Inverse problem 

In the linear approximation, the measured potential 
changes are assumed proportional to the conductivity 
changes. The images were reconstructed solving the 
normal matrix equation: 

A.Act = b with A = A T .A and b = A T Au (13) 
Act is the unknown vector of conductivity changes, 
Au is the vector of the measured potential changes and A 
the sensitivity matrix. The sensitivity matrix is ill- 
conditioned. The plot of singular values shows that the 
maximal rank of this matrix is 104 with 16 electrodes 
(Figure 5), which corresponds to the number of linearly 
independent measurements. Fan3 and fan4 show very 
similar sets of singular values. Experimental 
measurements confirmed that these two types of drive 
have equivalent performances so that fan4 was ignored 
in the following sections. The largest singular value of 
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adja, d=2, pu|| 2 =1.7 10 4 . N L =0.31% 



fan3. d=2. pu|| 2 =1.3 10 3 . N L =0 04% 





adja, d=3, ||Au|| 2 =2.0 10-\ N L =2.63% 



fan3, d=3, pu|| 2 =1.e 10 3 , N L =0.29% 



Figure 9 Images reconstructed using the same calculated data in Figure 8 with the addition of a constant gaussian 
noise voltage representing electrode and amplifier noise. The resulting noise level N L varied according 
to the sensitivity of each drive pattern and signal magnitude. The regularisation factor was X = 3. 16x1 0' 7 
for all images. 



adjacent drive is about four times smaller than that of 
fan3. 

Reconstruction method 

The equation was solved as an optimisation problem 
using Tikhonov's regularisation method, searching for a 
vector b minimising the functional F defined by: 



AArj • 



(14) 



with A = A A and b = A Au 

Symbol ||w|| p denote the p-norm of a vector w and X 
is the regularisation parameter. The optimal value of the 
regularisation parameter X was determined automatically 
for each data set using the "L-curve" procedure. The 
lower limit of X (10~ 15 ) was determined by successive 
trials using simulated data. For this value, the images 
reconstructed from simulated noiseless data could not be 
distinguished from reconstruction noise. The upper limit 
(10") corresponded to clearly excessive image 
smoothing producing lobes spreading over the entire 
image. In practice, the values found by the automatic In- 
curve procedure ranged roughly from 10" 9 to 10" 12 for 
noiseless simulated data and from 10" 6 to 10" 8 for 
experimental data. Regularised matrices were pre- 
calculated using three values of X per decade. Image 



reconstruction therefore consisted of matrix-vector 
products and calculation of the radius of curvature of the 
L-curve. The calculation was carried out using the 
circular mesh of Figure 2 using specific software written 
in Borland Delphi. The cartesian mesh for 3D plots of 
Figures 8 to 10 were obtained by linear interpolation. 
With a 1.6 MHz laptop PC computer the image was 
available in about 20 seconds. 



EXPERIMENTAL SETUP 

The experimental data set were collected using a 
bench system that comprised a 16-electrode mock-up 
probe, 50 mm in diameter, immersed in a tank filled with 
tap water modelling a uniform conductivity medium and 
the purpose-built experimental instrumentation [9]. The 
electrodes, made of brass, formed 50x2 mm 2 conducting 
stripes regularly spaced around the probe. The 
measurement frequency was 8 kHz [8]. This frequency 
was chosen for bench experiments as it ensured 
satisfactory compromise between the increase of 
electrode-medium interface impedance at low frequency 
and the onset of error due to stray capacitance with 
increasing measurement frequency. The magnitude of the 
measurement current was 1mA pp. The contact 
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adja, PVC rod, *=10 8 , (-) 



fan3, PVC rod, A = 1fJ- 8 , (-) 




adja, simulated, ^=10- 12 , (-) fan3, simulated, X=10 12 , (-) 



Figure 10 Images from experimental data obtained with a PVC cylindrical rod, 16 mm in diameter (normalised 
radius r = 0.32) and located at distance d = 3 (top row) compared to the images from simulated noiseless 
data of the same perturbation (bottom row). In this figure, the sign of the reconstructed values has been 
inversed to display the upward negative perturbations of conductivity. 



impedance was about 400 ohms per electrode in tap 
water. 

Tap water has been widely used in EIT as a uniform 
conductivity model, especially for feasibility studies and 
test of instrumentation, even though it does not have the 
same electric and dielectric properties as human body 
tissues. As a matter of fact, there is no really satisfactory 
model of the conductivity of cellular medium. The 
conductivity of tap water was 0.039 S/m in this study. 
Furthermore, the use of a liquid model makes it easier to 
place conductivity perturbations in the medium. 

The adaptable front-end system (Figure 6) 
comprised a controlled amplitude current source based 
on the circuit of Figure 3a, 16 input differential voltage 
amplifiers, 4 multiplexers for the 16 to 1 selection of 
electrodes. In this study, the demodulator produced a DC 
signal proportional to the magnitude of the measured 
voltage. This voltage was digitized with a resolution 
equivalent to 17 bits. This was achieved by means of a 
two-stage system reducing the DC offset of the signal 
using a D/A converter and sampling the amplified 
residual with a 12-bit AID converter. 



RESULTS 

Images reconstructed from noiseless calculated data 

Figure 8 shows the images reconstructed according 
to (14). The drive patterns fan3 and adjacent give similar 
images. For both types of drive, the resolution is better 
near the probe and deteriorates for increasing distance 
from the probe. 

Images reconstructed from calculated data with added 
noise 

Gaussian noise was added to the calculated data 
according to (11). The simulated real part of electrode 
contact impedance was 400 Q., simulating the value 
measured in vitro. The amplifier's input related noise 
figure was taken equal to 13 nV/VHz. Current noise was 
ignored. The resulting total noise figure was 
13.2nV/VHz. The bandwidth was B = 208. Figure 9 
shows the corresponding images. The conductivity 
perturbation does not distinguish from noise artefact, 
excepted for fan3 drive pattern at distance d = 2. This 
figure illustrates the influence of noise and the different 
susceptibility to noise of adjacent and fan3 drive patterns. 
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Figure 11 This animation displays successively the sensitivity maps of the 49 basic patterns from which any 
electrode pattern can be derived by symmetry and rotation. This animation does not display any time 
varying process. The purpose is to illustrate the influence of drive pattern on the sensitivity distribution 
around an EIE probe comprising 16 electrodes. The red colours show positive values of sensitivity and 
blue colours to negative ones. The highest magnitude values are near the electrode. The change is 10 dB 
from one colour level to the next. The central, background colour corresponds to low absolute values of 
sensitivity, either positive or negative. As there are 1 1 levels, the background corresponds to sensitivity 
smaller than 50 dB compared to the maximum. The maps were drawn varying the current injection pair 
of electrodes and keeping constant the sensing pair. The radius of the mapped zone is 3 times the radius 
of the probe. 



Images reconstructed from experimental data 

The experimental data sets were collected using the 
set-up described above. The perturbation used in this 
study was an insulating PVC, 16 mm in diameter 
(normalised radius of 0.32) and located at distance d = 3 
from the axis of the probe. Figure 10 shows the images 
reconstructed from the in vitro data and from calculated 
noiseless data modelling the same perturbation. 

DISCUSSION 

The above data and design parameters confirm that 
EIE is particularly sensitive to noise. This is mainly due 
to the rapid decrease of sensitivity with increasing 
distance from the probe due to the simultaneous 
reduction of the lead fields of current and voltage 
electrodes. The results obtained in this study are 
compatible with prostate size. The noise of the current 
injection circuit is negligible compared to the other 
sources of noise. The predominant source of noise is the 
input related noise of the voltage amplifier. The second 
largest source of noise is electrode noise. This gives 
particular importance to electrode technology in a 
miniaturised probe. 

For tissue characterisation, there is no particular 
need for high data acquisition speed. The parameter 



bandwidth (B) in (9) accounts for measurement time. 
The above numerical examples were based on the rate of 
1 frame per second with 208 values per frame. The 
increase of the total acquisition time by, for instance, a 
factor of 16, would be practically acceptable and would 
increase by 12dB the signal-to-noise ratio. With the limit 
fall-off in 1/d 4 of the sensitivity this would theoretically 
increase by two the sensitivity range, reduce the 
measurement current or improve the quality of the 
reconstructed images. One possible technique would be 
the averaging of 16 successive images. The use of a 
series of images would also enable the detection of 
transient artefact during data acquisition. 

The limit size for the detection of a tumour depends 
on the conductivity of the medium, the conductivity 
contrast of the tumour, the sensitivity distribution of the 
used drive pattern, the measurement noise and the 
magnitude of the applied current. There is therefore no 
unique answer to the question of limit size for detection. 
However, the plots in Figure 4 enable the calculation of 
estimates under the conditions of example numerical 
application. These plots correspond to the limit 
measurement conditions for the pixels of the mesh of 
Figure 2. Straightforward considerations based on 
equation 9 yield a value for given current magnitude and 
tumour location. Table 2 shows the values obtained with 
a conductivity ratio Ao/o = 10% for 1mA and 10mA 
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currents at the distance of three times the radius of the 
probe, assumed as typical range for EIE measurements. 
For extrapolation of the figures of Table 2 to 
conductivity changes than 0.1, the conductivity ratio 
q=Ao/o should be replaced with the quantity 2q/(q+2) 
ratio in equation 9 to account for non-linearity [35]. 

The two selected bipolar drive patterns were 
adjacent and fan3. Adjacent drive has widely been used 
in EIT due to its minimal number of measurements and 
its full ranked sensitivity matrix. Fan3 was selected in 
previous studies for its sensitivity range that was found 
larger than that of the other tested patterns [9,26]. The 
required number of measurements in fan3 is twice that of 
adjacent drive and the sensitivity matrix is not full rank. 
However, the additional linearly dependent data can be 
seen as an averaging improving signal-to-noise ratio in 
the same way as the full set of 208 measurements used in 
EIT to compensate for reciprocity error. The 
reconstructed images indicate that adjacent and fan3 give 
images of similar quality in absence of noise. Adjacent 
drive requires a four time larger measurement current for 
a given signal-to-noise ratio (Figure 4). This can be 
compensated by sufficient current magnitude and 
measurement time. Under these conditions, the above 
noise equation shows that both types of drive can be used. 
Preference would then be given to adjacent drive due to 
its better matrix conditioning. In practice, however, the 
use of small size electrodes could potentially limit the 
actual magnitude of the applied current below the 
maximal limit of 10 mA rms. Furthermore, the 
impedance of electrodes and the output swing of the 
current source can also limit the magnitude of the 
measurement current. Table 2 shows that fan3 enable the 
detection of larger tumour than fan3 for given 
measurement conditions. This argument may be decisive 
in selecting the drive pattern to be implemented. In any 
case, electrode technology will be crucial in the design of 
probes and that, finally, EIE seems more appropriate to 
tissue characterisation than to high speed imaging. 

CONCLUSION 

The simulation of operating conditions enabled the 
quantification of the magnitude of the measurement 
current ensuring appropriate signal-noise ratio. 
Measurement current of about 1 mA satisfying the noise 
condition derived in this study is feasible in practice. 
This study also showed that fan3 and adjacent drive 
patterns give similar results with noiseless data, but that 
adjacent drive requires significantly higher measurement 
current than fan3. The equation derived in this study 
enables the specification of the hardware system given 
the operating condition of a given application. 
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Table 1 Limit values of the 2D -calculated and experimental geometry factors for fan3 and adjacent drive 
patterns 





fan3 


fan3 


adjacent 


adjacent 




calculated 


measured 


calculated 


measured 


§max 


0.333 


0.254 


0.0958 


0.151 




0.125 


0.123 


0.0124 


0.00945 



Table 2 Limit diameter of a conductivity perturbation (Ao/o = 10%) located at three times the radius of a probe 
under the example measurement conditions. Figures in percent are relative to the diameter of the probe. 
Figures in mm are for a probe 7 mm in diameter 





adjacent 


fan3 


adjacent O = 7 mm 


fan3 O = 7 mm 


1 mA 


25% 


8.4% 


1.8 mm 


0.6 mm 


10 mA 


7.9% 


2.7% 


0.6 mm 


0.2 mm 



